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Abstract 

Apparent regularities in fermion masses and mixings are often associated with physics at 
a high flavor scale, especially in the context of discrete flavor symmetries. One of the main 
reasons for that is that the correct vacuum alignment requires usually some high scale 
mechanism to be phenomenologically acceptable. Contrary to this expectation, we present 
in this paper a renormalizable radiative neutrino mass model with an A^ flavor symmetry 
in the lepton sector, which is broken at the electroweak scale. For that we use a novel way 
to achieve the VEV alignment via an extended symmetry in the flavon potential proposed 
before by two of the authors. We discuss various phenomenological consequences for the 
lepton sector and show how the remnants of the flavor symmetry suppress large lepton 
flavor violating processes. The model naturally includes a dark matter candidate, whose 
phenomenology we outline. Finally, we sketch possible extensions to the quark sector and 
discuss its implications for the LHC, especially how an enhanced diphoton rate for the 
resonance at 125 GeV can be explained within this model. 
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1 Introduction 



Much progress has been achieved in the field of particle physics during the last year. First, 
the last missing mixing angle #13 of the Standard Model(SM) with massive neutrinos has 
been measured [ ] to be 8° after first hints in 2011 [2] and recently an excess consistent with 
a SM Higgs has been observed at 126.0 ± 0.4(stat.) ± 0.4(sys.) GeV by ATLAS [ ; ] and at 
125.8 ± 0.4(stat.) ± 0.4(sys.) GeV by CMS [5; 6]. 

Let us first discuss the implications of the large mixing angle #13. Much of the work 
in the neutrino sector has been aimed at explaining tiny values of #13 as deviations from 
a tri-bi-maximal(TBM) mixing structure using flavor symmetries, but this scenario is now 
implausible due to the sizable value of #13 . So maybe there exists no symmetry that is 
connected to the regularities in the fermion parameters. It could rather be that mixing angles 
are determined at a high scale from some (quasi-)random mechanism. Indeed if one randomly 
draws unitary 3x3 matrices with a probability measure given by the Haar measure of U(3), 
i.e. the unique measure that is invariant under a change of basis for the three generations, 
one finds a probability of 44% for nature to have taken a more "unusual" choice [9]. This 
cannot be interpreted, however, as an indication in favor of anarchy [10; 11], as the sample (3 
mixing angles and one mass ratio) is clearly too small to reconstruct the probability measure 
to any degree of certainty [ 2]. The only statement one can make is that the (very limited) 
data cannot rule out the anarchy hypothesis. For any values of the mixing angles one can 
always find a flavor model that is in better agreement with the data 2 . 

Another option is that flavor symmetries are realized in a different way. One route is to 
think of solutions that do not predict TBM. Such models are usually implemented at high 
scales and give precise predictions for the leptonic mixing angles in the experimentally allowed 
regions. These models might then be falsified in the same way that the models that give tri- 
bi- maximal mixing have been ruled out, i.e. by a further refinement of the experimental 
determination of these angles. This seems to be the only fruitful direction for models that 
explain flavor at high energy scales such as the see-saw or GUT scale, because mixing angles 
are generically the only experimentally testable predictions of such models. An example are 
models based on A(96) [14-16]. 

An interesting question is if flavor symmetries could even be realized at low scales [17-28]. 
If this is viable then such models can be tested by additional observables. Such observables 
typically include rare lepton flavor violating(LFV) decays of leptons and mesons and -ideally- 
a direct experimental access to the very fields that mediate the flavor symmetry breaking. 
Typically such models will feature extended Higgs sectors but will not uniquely determine the 
mixing angles; they will, however, rather give relations among the deviations from patterns 
such as tri-bi-maximal mixing. 

In this work we implement a model based on a flavor symmetry at the electroweak scale 
and show that it leads to additional phenomenological effects that will become testable in the 
future. 

Many authors have noted that quite generic corrections in the neutrino sector within 

1 See also the recent global fits to neutrino oscillations [7; 8]. 

2 This can be done without increasing the degree of complexity of the model [13]. 
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models based on the symmetry group A4 [17; 2 1] may lead to corrections to the leptonic 
mixing angles that are in agreement with the experimental data [30]. It should be noted, 
however, that all such models need to break the symmetry group in a particular way to 
different subgroups in the charged lepton and neutrino sectors, a vacuum configuration that 
cannot be obtained from a straightforward minimization of the potential but something that 
rather needs a special dynamical mechanism to achieve it. The two most commonly used 
mechanisms are either based on (i) continuous R-symmetries in supersymmetry or (ii) extra- 
dimensions. The supersymmetric models only work in the limit where the supersymmetry 
breaking scale is below the scale of flavor symmetry breaking and the scale of flavor symmetry 
breaking is thus unobservabely high. The only experimentally verifiable prediction of such 
models seem to be correlations in the deviations from TBM, e.g. the trimaximal mixing 
pattern that predicts a ~ — ^rcos<5 and s ~ where [31] 

r 1 1 

sin6»i3 = -^=, sin(9i2 = -j= (1 + s) , sin0 2 3 = ^= (1 + a) . (1.1) 

As these models tend to be quite baroque, since they involve a large number of driving 
fields etc., this one prediction seems to be a rather poor showing compared to the model 
building effort involved. The second possibility using extra-dimensions is possible, but we 
here want to focus on non-supersymmetric models in four dimension, a possibility that seems 
to be favored by the experimental data. 

In an earlier work [32], two of us(MH and MS) have studied a possibility to obtain the 
vacuum alignment that seems to be quite close to the general spirit of discrete flavor model 
building: to get a predictive model one needs to have unbroken discrete remnant symmetries 
in the neutrino and charged lepton sectors that do not commute with each other. The way this 
is usually realized is that symmetries, particle content and vacuum structure are selected such 
that these remnant symmetries are realized as accidental symmetries of the leading-order (LO) 
mass matrices that are broken at next-to-leading order. The idea of [32], following earlier work 
by Babu and Gabriel [33] , was to have an extended flavor group G and a particle content such 
that an accidental symmetry G x A4 arises at the renormalizable level in the scalar potential 
that allows for the desired vacuum configuration. More precisely, the scalars \% breaking the 
extended flavor group G in the charged lepton sector transform under the A4 only, while 
the scalars fa breaking the extended flavor group G transform under the full group G. The 
symmetry is chosen such that the scalar potential does not contain operators with non-trivial 
contractions of Xi with (pi, i.e. there are only contractions of the form (XiXi)l 1 • The 

smallest symmetry group that realizes such a structure for A<± is G = Qs x A4 and since the 
symmetry breaking does not need any special additional ingredient, there is no immediate 
theoretical obstacle to have the flavor symmetry breaking scale at the electroweak scale. 

In this work we therefore implement the model in [32] at the electroweak scale. The 
outline of the paper is as follows: We introduce the model and the exact symmetry breaking 
pattern in sec. 2. Without the introduction of any additional symmetries, we show in sec. 3 
that neutrino masses are generated at one-loop level and that the neutrino mass matrix is 
determined by five physical parameters. We discuss the predictions for neutrino oscillation 
observables that follow from this structure. In sec. 4, we discuss constraints from lepton-flavor 
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violating rare decays. In sec. 5, we show that the model contains a dark matter candidate and 
discuss its phenomenology. Three simple possible extensions to the quark sector are discussed 
in sec. 6. Finally, in sec. 7 we discuss the implications of direct collider searches and the the 
recent observation of a Higgs-like boson at 125 GeV and then we conclude in sec. 8. 



2 Model and Symmetry Breaking 

We utilize the symmetry Qg x A4 introduced in [32], which allows for natural vacuum align- 
ment, and implement a model describing the lepton sector at the electroweak scale. Hence, we 
promote the flavon fields of [32] that couple to the charged lepton sector to EW Higgs doublets. 
The particle content of the lepton sector is given in Tab. 1. The vacuum configuration 

(Xi) = ( ^ 1 , <0i) = ^=(a, a, b, -b) T , (0 2 ) = -t=(c, c, d, -d) T (2.1) 

can be naturally obtained from the most general scalar potential following the discussion in 
[32]. As the discussion is very similar to the one given there, we relegate it to appendix A.l, 
where also the scalar mass spectrum is discussed. However, let us briefly recall the salient 
features of the vacuum expectation value (VEV) configuration (2.1): the scalar singlets 4>\ and 
4>2 break the symmetry group to the subgroup (5|5 2 = = Z2 and the EW doublets x break 
the discrete symmetry group down to the subgroup (T\T S = E^ = Z3, while simultaneously 
breaking the electroweak gauge group SU(2)i x U(l)y down to the electromagnetic U(l) em - 
The normalization is chosen such that J^vf = v 2 = {y/2GF) 1 = (246 GeV) 2 , in accordance 
with our earlier definition. Because of the unbroken Z3 symmetry in the charged lepton 
sector, it is useful to go to a basis [17; 22-24] 

(ff,^>") T = 4x~(l,wV), (L e ,L^,L T ) T =tf T L~ (l,w 2 ,u;), (2.2) 

where this symmetry is represented diagonally and Cl? is defined by 
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(2.3) 



We have indicated the transformation properties under the unbroken subgroup (T) = Z3 
under which (e c ,/i c ,r c ) transform as (l,w,u; 2 ). This has been denoted flavor triality in [22]. 
In this basis the vacuum configuration (2.1) implies that only the field H acquires a VEV 
(H) = (0, v/y/2) , while ip' and <p" are inert doublets (and thus do not obtain a VEV). The 
potential for the electroweak doublets \ is given by 

V x (x) = (4x j X + Yl V(x t x)r(x t X)r* + A x Alm [(x+x)^ (x t x)3 1A l . ( 2 - 4 ) 
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and after symmetry breaking the nine physical scalars contained in \ arrange themselves in 
the following multiplets under the remnant U(l) x Z3 symmetry. There is one real scalar 
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Table 1: Particle content of the minimal model that realizes flavor symmetry breaking at the electroweak 
scale. The flavon \ contains the Higgs field and ties the electroweak to the flavor breaking scale. The scalars 
r]i and fermionic multiplet S are needed for one-loop generation of neutrino masses. 



3A xll + ^A x3l , s K (2.5) 



h = V2ReH° with mass 

2 2 
h 9 

that plays the role of the Standard Model Higgs. Note that since this scalar is a complete 
singlet under all remnant symmetries, it can in principle mix with components of 0i and 02 
that transform in the same way. This is discussed in detail in Eq. (A. 9) in the appendix and 
in the following we will for the most part assume the mixing to be small enough to treat h as 
a mass eigenstate. 

The next four degrees of freedom are in the charged scalars (p' + and ip" + that transform 
as (l,w 2 ) and (l,w) under U(l) x Z$, respectively, and have the masses 

2 2 

m l'+ = ^ (-2v / 3A x 3l , s - A xA ) , "V+ = (" 2 ^ A x 3 liS + A x a) • (2.6) 

The final four real scalars sit in the two complex neutral scalars tp'° and tp"° , that both 
transform as (0,u; 2 ) and the mass eigenstates are given by the neutral scalars 

= ( cosa sinQ 1 1 V I (2.7) 

\ <f J \ — sin a cos a ' 1 ' 
the mixing angle a and their masses may be written as 
, 6A x 12 + ^ (3A X 3 A + A x 3 s ) 

tan 2a = ■ , (2. 

6A xA ; 

,2 1 2 o+„„^o„\ /™2 2 \ 1 2 1 2 u A X 3i,a 
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m 01 + m^ 2 = -2 tan(2o) ( - m v ,+ ) + + m^, + -j=^- , (2.9) 



| 2 = 2 |sec(2a)| ui\, l+ - m",+ . (2.10) 



The mass spectra for the other scalars can be found in appendix A. 2. Two comments are 
in order here: (i) in the potential (2.4) there is only one mass term for the three doublets. 
Using the minimization conditions, the mass term can be swapped for the Higgs VEV v and 
therefore (ii) all of the squared scalar masses are given as a product of dimensionless scalar 
couplings times v 2 . The additional scalar masses may therefore not be arbitrarily large. Note 



5 



that in usual multi-Higgs doublet models each doublet has its own mass term and therefore 
there is always a decoupling limit where all non-SM particles are unobservabely heavy. Such 
a setup is therefore directly testable at colliders, as we will study in sec. 7. However, before 
discussing this, we show that the model accomplishes (i) the description of the (lepton) flavor 
structure in terms of a small number of parameters and (ii) the protection against bounds on 
new physics from flavor observables such as lepton flavor violating processes. 

3 Lepton Flavor Structure 

In this section we discuss the one-loop generation of neutrino masses and phenomenological 
implications of the predicted flavor structure. 

3.1 Lepton Masses 

The charged lepton sector is described by 

- C e = y e Lxe c + y^Lxfi + y T Lxr c + h.c. , (3.1) 

where x = l(T 2X* an d here and in the following we do not specifically indicate the contractions 
if there is only one invariant that can be formed out of the particle content of the operator. 
In the physical basis of Eq. (2.2) this term reads 

-C e =H (y e L e e c + y^L^fi + y T L T r c ) + {y e L^e c + y^L T ^ c + y T L e r c ) 

+ <p" (y e L T e c + y ll L e fi c + y T L )l T c ) + h.c. (3.2) 

and we thus see that H couples diagonally to leptons while cp' and if" do not. Note that here 
the mass terms are generated by dimension four Yukawa couplings and there is therefore no 
need for a complicated UV completion. The mass matrix is thus given by 

v 

M E = -j=£l* T &\a.g{y e , y M , y T ) , (3.3) 

with £1? given in Eq. (2.3). Neutrino masses are generated at one loop level, through the 
interactions with the fermionic singlets S and the scalar doublets 77, as shown in Fig. 1. The 
couplings of S are given by 

Cu = hiLrnS + h 2 Lrj 2 S + VZM S SS + h.c. . (3.4) 

The factor of \/3 cancels a factor coming from the normalization of Clebsch- Gordon coeffi- 
cients. In order to calculate the neutrino mass matrix, we have to determine the mass matrix 
of the neutral components of 771, i] 2 and 773. To shorten the notation we define the doublet fjj 
to be the J-th component of the 9 component vector f) = (7/1,772,773) and real scalar field 77° 
to be the A;-th component of (v^Re?) ' : ^/2\mfp). Besides the direct mass terms 

f)2T/( 2 ) 

K)« = lw with ^ a) =£^M*. ( 3 - 5 ) 

J ''<■ i=l,2,3 
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Figure 1: Neutrino mass generation at one loop. 



there are couplings which give off-diagonal contributions 




to the mass matrix. Such interactions are needed to generate neutrino masses and the relevant 
ones can be determined from symmetry considerations 3 . Any contribution to neutrino mass 
has to be proportional to 

• Ms, which breaks the generalized lepton number L — > e ia L, S —> e~ ia S 

• either of the couplings Ai or A2, defined by 4 

v v,x = A i(x T c r 2'?x)l 1 (^fc r 2'??/3)i 1 + A 2 e lcVA (x T '720 ; x)3 1 ( 7 ?2 ro '20 ; ??3)3 1 + h.c. , (3.7) 

which break the generalized lepton number L — > e ia L,rji — > e~ ia r]i, 

• and A3 or A4 defined by 5 

Vr,j =\ 3 (4>i4>2)l 1 + A 4 (0102)3 1 (4%)^ + h - c - » ( 3 - 8 ) 

which couples to the ^-breaking VEV of (fa. 

The built-in multiple protection of the neutrino mass operator thus necessitates the large 
number of couplings involved in neutrino mass generation, and thus a large potential for 
suppression beyond the naive factor of l/(167r 2 ) from the loop integral. For simplicity, we 

3 The complete expression for SVrjf* can be found in the appendix in Eq. (A. 10). Here, we only present the 
parts which are relevant for neutrino masses. 

4 We can set a number of complex parameters real by phase redefinitions. We set y e , y M , y T , hi, /12, Ms, 

Ai, A3, A4 real by rotating £ c , L, r)2, S, x, ill, T)3, respectively, and display the phase of A2 explicitly. 

5 The contractions {x T cr 2^X)\ vanish in the vacuum given in Eq. (2.1) and thus do not contribute to 

2,3 

the masses here, because the Z3 symmetry generated by T is conserved by {%)■ 
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assume that the direct mass terms M{ dominate over all other contributions; this is in fact 
a necessary condition to have a predictive theory of flavor. Hence, we can approximate the 
propagator as 



k 2 - (M 2 + 5M 2 )] = (k 2 - M^o)" 1 + (k 2 - M^yHM^k 2 - M 2 )~\ (3.9) 



where M 2 is diagonal, and treat the mixing between the different components of rji by mass 
insertions 5M 2 . The evaluation of the one loop diagram leads to 



3 18 

(M„U = Y1 E h aiI h^jM s (dM 2 ) iM {SM 2 
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,M S ) (3.10) 



where the Yukawa couplings h^j depend on the two couplings h\ 2 given in Eq. (3.4) via 



h-akj 



dL a dS k df]j 



(3.11) 



and the loop integral is given by 6 

I(m 1 ,m 2 ,m 3 ,m i ) 



JL V— — 
6^ Z^fTV. 



im 2 log ( ^ 



(3.12) 



16n' 2 ^ n fc ^ (mi 2 - m k 2 ) 
Evaluation of the sums leads to the following flavor structure of the neutrino mass matrix: 



Mu 



V 



a + be iax d + ee iax 
a 



(3.13) 



where the four real coefficients are given by 

1 
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h\ A 3 \iv 2 (ac + bd) M s I (Mi, Mi, M 3 , M s ) , 
h x h 2 A 4 Xiv 2 {be - ad) M s I(M X ,M 2 ,M 3 , M s ), 
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h\ A 4 X 2 v 2 (be - ad) M s I(M 2 , M 2 ,M 3 ,M S ), 
h x h 2 X 3 X 2 v 2 (ac + bd) M s I(M X , M 2 , M 3 , M s ). 



(3.14a) 

(3.14b) 

(3.14c) 
(3.14d) 



Hence, neutrino masses are suppressed by one insertion of the EW breaking VEV Ai (x 2 ) /Mq, 
with Mq being the largest mass of the particles in the loop Mo ~ maXi=i,2,3,S Mi, and one mass 
insertion of the flavor breaking VEV A2 (4>i4> 2 ) /Mq. A phenomenologically viable neutrino 



6 Note that the renormalization scale fi drops out of the sum; it is displayed here to make the symmetric 
structure of the expression explicit, while keeping the argument of the logarithm dimensionless. 
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mass scale is obtained for e.g. M ~ O(TeV), (x) , (fc) ~ 0(100 GeV) and h h A* ~ 0(0.01 - 
0.1). The next-to-leading order corrections are suppressed by Ai (x 2 ) /Mq or A2 {4>i4>2} /Mq, 
which amounts to an 0(0.0001 — 0.001) correction for our typical values and can be neglected 
to a good approximation. 

The neutrino mass elements correspond to the following operators: 

• a : (L T a 2 aL)^ (x T T2Tx)i (^1^2)^ 

• d : (L T a 2 aL) 3 ^ (x T <72<?x) 1 (01^2)3^ 

• e : (L T a 2 aL)* ix T ^x)^ (^l^h., 

• Ei w<_1 {L T a 2 aL) lM 

The fact that only the combination shown in the last line contributes to neutrino masses is 
due to the UV completion presented here. In a general theory one might have all operators 
present, thereby reducing the predictability of the theory. 



3.2 Phenomenological Implications 



As the neutrino mass matrix is described by five physical real parameters, there are four 
predictions in the lepton sector at leading order. They can easily be read off from Eq. (3.13) 
in terms of matrix elements, but the expressions in terms of mixing parameters are non- 
trivial. In the flavor basis, where the charged lepton mass matrix is diagonal, the neutrino 
mass matrix is given by 



Mj l 
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(3.15) 
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and it is instructive to look at the neutrino mass matrix in the tri-bimaximal basis M, 
U]j PS MJ 1 Uhps, i.e., 
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(3.16) 



We will first discuss limiting cases analytically and then perform a numerical analysis of the 
general neutrino mass matrix. In the limit of | A2 1^ 2 — > 0, both b and e vanish and we obtain 
tri-bimaximal mixing 
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Figure 2: The deviations from tri-bimaximal mixing of the form U = UhpsVxi and U = UhpsUis generated 
by the angle 9 defined in Eq. (3.18). The yellow point represents TBM, the continuous lines give the deviations 
from TBM with the angle 9 given by the color codes in the top right corner for 5 = ^ ^ for n = 0, . . . , 5, where 
n — is the outermost parabola etc. The one, two and three sigma regions of a recent global fit [7] are 
indicated by dotted, dashed and continuous contours, respectively. 



From Eq. (3.16) we can read off that switching on e / while keeping 6 = results in a 
correction to the PMNS matrix of the form U = UhpsUi2{9i2)P with ?7i2(#i2) denoting the 
unitary matrix 

/ ci2 -s 12 e- iSl2 \ 
Ui 2 (6) = s 12 e iS ^ c 12 , (3.18) 



with C12 = cos #i2, s\2 = sin #12 and P being an arbitrary phase matrix. In the standard 
parameterization of the PMNS matrix [34] with the 1-2 rotation to the right, this 1-2 correction 
only affects the solar angle, while maintaining the predictions of a maximal atmospheric 
and vanishing reactor angle. Since large corrections to this angle are not allowed, in the 
phenomenologically acceptable parameter space the relations e <^ b,a,d should hold. 

On the other hand, if we take b ^ while e = 0, we see from Eq. (3.16) that this requires a 
1-3 correction U = U h psUi^{9i^)P , where U\^{6i^), analogous to C/i2(#i2), denotes a complex 
rotation in the 1-3 plane. This correction is of the trimaximal mixing [30; 31; 35-39] form, 
which can perturb TBM back into agreement with experiment. The effect of the various 
deviations from TBM is illustrated in Fig. 2. 

To gain an analytical understanding of how the additional parameters affect the mixing 
angles, we can perform a perturbative analysis in the limit of small e and therefore small 
| sin 2 #12 — ||. The PMNS matrix can be described by UhpsUiz{Qi3)Ui2{0i2)P , where #12 and 
#13 are small in the phenomenologically interesting region and the Majorana phases are given 
by P = diag(e iai / 2 , 1, e ia ^ 2 ). Hence, we can permute the matrices ?7i2 and C/13 and we define 
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r\i = sin#ij cos<5ij and tu = tan<5ij, which evaluate to 
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(3.19a) 
(3.19b) 



where Am^ is the leading order solar mass squared difference, i.e. neglecting the small 
corrections of 7*13 and e. The phases of the matrix P are given by 



tan ol\ 
tan 012 



26ri3 cos a\ — sinct\(2&ri3ii3 + b + 2e) 
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2(d + 6ri3sinaA — d) + cosaA(— 26ri3ii3 + 6 — 2e) 



(3.20a) 
(3.20b) 



Similar to [31], we can parameterize the leptonic mixing matrix in terms of deviations from 
the tri-bimaximal mixing angles as defined in Eq. (1.1). The Dirac CP phase 5c p is undefined 
in the tri-bimaximal mixing limit and we leave it free and do not expand in it. Besides the 
contributions of ot\ 3 to the Majorana phases (p± 2 in the standard parameterization, there are 
also small corrections 6ipi2 from the matrices ^12(^12) and ^13(^13) resulting in 



<Pi = «i - "3 + Scpi, 



and 



ip 2 = ir — as + 5cf2 



(3.21) 



This expansion leads to the following form of the PMNS matrix 

f s+i8tp\—2 
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(3.22) 



Equating the expanded form of U to UhpsU\z(Qiz)Ui2(Qi2)P determines all free parameters 
s,r,a,S,S(fi,S(p2 as well as some corrections to unphysical phases, which we suppressed for 
simplicity. The first order deviations from the mixing angles are 



s = — v2ri 2 £i2, rcos<5 = - 

and the CP phases are given by 

tan Sep = tan ^13 ipi = a\ — as — 2\/2ri2 



2r i3 _ ri3 

x/3' "V3' 



yj 2 = 7r - a 3 - 2\/2ri2 . 



(3.23) 



(3.24) 



Following [31], we can derive a sum rule, which relates the deviations of the atmospheric 
mixing angle with the ones of the reactor mixing angle 



1 



r cos Sep ■ 



(3.25) 



The masses are determined by 
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Figure 3: Dependence of the reactor angle #13 on the atmospheric mixing angle 823- The various color 
codings are given next to each scatter plot. Top left: For sin 2 623 < l/2( sin 2 823 > 1/2 ) the model predicts 
Sep = 0,2tv(8cp = 7r). Top right: The scatterplot shows a band structure in sin 2 #12. Bottom left: For the 
points in the experimentally allowed region, b has to be of similar size as a, d. Bottom right: For the points in 
the experimentally allowed region, e has to be of approximately one order of magnitude smaller than a, d. The 
one, two and three sigma regions of Ref. [7] are again indicated by dotted, dashed and continuous contours, 
respectively. 



. b 2 

m\ = a 2 + b(d + d) cos a\ + 2dd + — + a 2 , m\ = a 2 , (3.26) 

h 2 

m 3 = a + b (a — d) cos a\ — 2dd + — + a , 

to leading order in the small mixings ri3, ri2, and the leading order ratio of mass squared 
differences is given by 



Amli 4a (2d + b cos a\) + 4d(d + b cos at\) + b 2 
Am 2 i2 4 a (2d - b cos a A ) - 4d (d - b cos a A ) - b 2 

At next-to leading order, rrt\ and 771,3 receive corrections 



(3.27) 



Sm\ = b(2ris(d + d) sinaA + bri S ti 3 + e) + 2(d + d) cosQA(6ri 3 ti 3 + e) , (3.28a) 
(5m| = — 6(— 2ri3(a — d) sina^ + br 13 ti 3 + e) — 2(a — d) cos 0:^(^13*13 + e) . (3.28b) 

To illustrate our findings numerically, we have performed a numerical scan over the model's 
parameter space. We have randomly drawn values for the model parameters of order unity, 
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Figure 4: Numerical evaluation of the approximate 
atmospheric sum rule (3.25). The numerical evalua- 
tion shows that the sum rule holds to a good degree 
of approximation. 
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Figure 5: Expectation for the effective mass of neu- 
trinolcss double beta decay. The pink points lie 
within the 3 sigma region for all oscillation para- 
meters. The points with color coding lie within the 3 
sigma range for all observables except #13. 



assuming a Gaussian distribution with an expectation value of one and a variance of 0.5. The 
plots in Fig. 3 show the relation between the atmospheric mixing angle 623 and the reactor 
angle #13 . From the bottom two plots one can read off that b is of the same order as a and d for 
the experimentally measured #13 while e has to be about one order of magnitude smaller. The 
color codings of the two top panels show the mixing parameters 5c p and sin 2 #12. Clearly, 
the model is predictive: if sin 2 #23 is found to be close to the best fit point in the octant 
with sin 2 #23 < 1/2, the prediction for the CP phase is 5cp = 0, 2tt while for sin 2 #23 > 1/2 
it is predicted to be 5c p = vr. To establish the correlation with sin 2 9 12 shown in the top 
right panel, a precision determination of all the mixing angles is needed. In Fig. 4, as a 
consistency check of our analytical expressions, the atmospheric sum rule (3.25) is shown 
for the points obtained in the numerical scan. The color coding gives an indication of the 
magnitude of deviations from TBM and for small values the approximate relation is fulfilled 
to good accuracy. 

Finally, let us comment on the predictions for neutrinoless double beta decay. As can be 
read off from Eq. (3.15), the effective Majorana mass of the electron neutrino is given by 



2d 



+ 2e + 



(3.29) 



which can be expressed in terms of physical parameters as 



\m e 



mi 



2mi — 777-2 



2mi + 2?7l2 2<5(£>1 777,1 — 6<p2lTl2 

s — i 

2m\ — 777.2 2777,1 — 7772 



(3.30) 



As the additional neutral fermions S do not mix with neutrinos, there is no additional con- 
tribution due to the heavy singlet, like in Ma's scotogenic model [18; 19]. In Fig. 5 we show 



13 



the predicted range for the effective Majorana mass of the electron neutrino. As can be seen, 
the scan of parameters prefers moderately large values of the absolute mass scale, however, 
the effective Majorana mass of the electron neutrino can become small or even vanish. 



4 Lepton Flavor Violation 

In models with radiative neutrino mass generation, generally the particles in the loop can also 
mediate flavor changing processes, in particular lepton flavor violating rare decays. Before we 
enter into a detailed discussion of the various processes, we want to remind the reader about 
the remnant Z3 symmetry in the charged lepton sector 

(#,V?V) ~ (l,w 2 ,w), (L e ,L M ,L T )~(l,u; 2 ,u;), (e c , //, r c ) ~ (1, u, u 2 ), (4.1) 

which suppresses several LFV rare decays. If the remnant Z3 would be a symmetry of the 
whole Lagrangian, only the following LFV rare decays 

t + — > n + fi + e~ and r + — > e + e + /i~ 

and their charged conjugates would be allowed. All other decays can only proceed through 
a coupling to the Z3 breaking VEVs of the neutrino sector. Those decays are naturally 
suppressed and the symmetry thus protects the model from large constraints. At first, we 
will discuss the radiative LFV rare decays — > lj"f in sec. 4.1, focusing on the experimentally 
most well studied process, namely the process \x — > ej. In sec. 4.2, we discuss the LFV rare 
decays with purely leptonic final states, which are allowed at tree level, but suppressed by a 
three-body final state. Finally, we calculate the anomalous magnetic moment of the muon 
and compare it to experiment in sec. 4.3. 

4.1 Radiative LFV Decays /j — > ^7 

Let us first discuss the process of type U — > Ijj using an effective field theory approach. Such 
processes are described by effective operators of the form [40; 41] 

La^F^£ c H/M 2 ~ (3 1; 1) , (4.2) 

which transforms in the same way as the mass term under the flavor symmetry. It thus has 
to be multiplied by flavons to form an invariant. As we already mentioned, the remnant 
Z% symmetry in the charged lepton sector forbids all radiative LFV rare decays. Hence, the 
effective operator in Eq. (4.2) has to involve VEVs of the neutrino sector in order to lead to 
non-vanishing decay rates. The lowest order operators that can multiply the mentioned LFV 
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operator in the flavor basis read 



] / ^t)s 1 ) = l(ab(b 2 -a 2 ))(l,l,lf, (4.3a) 

4 ( ) = l ~ c ')) C 1 . ^ ( 43b ) 

4 (fafrl)^ ) = ^ ( a& ( c2 - d2 )) (i, !> 1 ) T ) ( 4 - 3c ) 

4 ( (^1^-1)3 ) = l {cd(a 2 - b 2 )) (1,1, 1) T . (4.3d) 



There can be more than one contraction, but in the vacuum they all result in these expressions. 
The lowest order effective operators thus all give contributions that can be written as 



C eff = i^t cT H^a^F^ML + h.c. with M 



( a\ a\ a\ \ ,, AS 

«2 «2 «2 ] ^ (U) 
y a 3 a 3 a 3 



where on are dimensionless couplings that should (naturally) be of order one and the mass 
scale M is the suppression scale of the higher dimensional operators. Note that the structure 
of flavor symmetry breaking in the neutrino sector is encoded in Ai. The symmetry thus 
automatically leads to a large suppression. From this matrix the LFV transition amplitudes 
can be determined as [40] 

and the magnetic dipole moments etj and electric dipole moments of the charged leptons 
are given by [40] 

v v 

a; = 2rrii—= KeMu, di = e—= ImAia. (4.6) 

V2M 2 V2M 2 V ' 

Note that the matrix M has additional dominant contributions to the diagonal entries stem- 
ming from operators that involve \ instead of (^) 4 . Using only the observables \i — > ej, 
t — > [AJ and r — t- e7 as well as charged lepton electric and dipole moments, it is therefore 
very difficult to test the underlying symmetry pattern, but it can give important indications 
distinguishing different models. For example in this model one would expect - barring the 
possibility of fine-tuned cancellations among the a, - similar branching ratios for the LFV 
decays \i — > ej, r — > fi'y and r — > ej, as was also found in SUSY A± models [40; 42]. 

In the following, we will focus on jj, — > &y, which is the most tightly constrained LFV rare 
decay. The leading contribution to /i — > e^f is given by the diagram depicted in Fig. 6a. It is 
similar to the neutrino mass diagram Fig. 1 in the last section. LFV rare decays mediated 
by the flavor violating EW doublets f'^ are suppressed by one more loop order because of 
the necessity to couple to the neutrino sector VEVs. Hence, they only show up at two loop 
order, as shown in Fig. 6b. We will therefore not consider this diagram further. 
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(a) Lowest order ;U — > e7 pro- 
cess mediated by 77's. 



(0102> '/ <01<fe) 



— ■» a » » ^ »^. e 

t 7 ( r ^^^^ 

(b) Lowest order fi — > e7 pro- 
cess mediated by </j'. 



Figure 6: Lowest order /i — > ej processes mediated by rj (left) and cp'^' (right). There has to a coupling to 
the VEVs (0i02) of the neutrino sector, which suppresses the amplitudes. 



Without any mass insertion along the r\ line, a one-loop diagram of this type evaluates 

to [19; 43; 44] 



Br(/i — > erf 



3 a 



64vr(G F m2)2 



(4.7) 



where m 2 Q = \ (Mf + Mf + Mf) and, using xj = {M 2 + )jj/ml and h akJ = gzM^Mj ■ 



C 



3 9 



jj 



=1 j=i 



and F 2 (t) 



1 -6t + 3t 2 + 2t 3 - 6t 2 lnt 
6(1 - ' 



In our model, we have C 2 = for the symmetry reasons given above and there have to be 
mass insertions to generate flavor violating interactions. Note that this is a welcome feature 
since LFV processes of this type severely constrain models that generate neutrino masses 
radiatively [. ]. This can be seen as the experimental constraint Br(/x — > erf < 2.4 -10 -12 [45] 
requires C 4 ~ 1.5 • 10 -8 for M5 = mo = 100 GeV. The flavor symmetry automatically reduces 



C by a factor 



SM 2 
v 



explicitly and we find 



In the limit 



Br(// — > erf 



SM 2 
v_ 

M 2 , 

7)1 



<C 1, the diagram 6a can be computed 



o 



16vr(G F m 2 ) 2 



(4- 



where 



C 2 



1 



m. 



E E " 



5M 2 



SM. 



KL 



h* eiL F A (M s ,Mj,M K ,M L ) 



(4.9) 



i=l J,K,L=1 

and -F4 is a dimensionless loop integral, which we only give in the limit of degenerate ry masses 



G 2 (t) = F 4 (M S = tmo, Mj = m , M K = m , M L = m ) 

= 4g ^ 2 \ ^ 12 [1 - 12t 2 - 36t 4 + 44t 6 + 3t 8 - 24(2t 2 + 3)t 4 Int] . 



(4.10) 
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The dimensionless functions F2 and G2 are plotted in Fig. 7. The explicit form of the sum 
in the expression (4.9) for C 2 is quite involved and will not be shown here, but it can be 
easily obtained using Eq. (A. 10) from the appendix. Here, we only comment on the generic 
size of the branching ratio. In general, the processes — > cy and the radiative neutrino mass 
diagram break different approximate symmetries and it is therefore not necessarily the case 
that the smallness of neutrino masses implies a small branching ratio. This is also the case 
here. For example from Eq. (3.14), one can read off that the smallness of neutrino mass could 
be due to very small values for Ai ~ A2 ~ 10 -9 , with all other couplings being order one. 
Then the dominant contributions to C 2 would be of the type 

G 2 (t) 



C z D 



1 



ml 432 



/i2A4(6c — ad) [— h\\%(ac + bd) + uj 2 h2\n{bc — ad)] 



(4.11) 



where we have again used the limit of degenerate 
and could in principle be of order one. However, if we 
stick to the parts of parameter space where the small- 
ness of neutrino mass is due to many moderately small 
couplings hi,\i ~ 0(0.01 - 0.1) and m ~ O(TeV), 
(x) , ~ O(100 GeV) (as discussed below (3.14)) in- 
stead of one very small coupling, the branching ratio 
is heavily suppressed by C A ~ (10 -9 — 10~ 13 ) 2 . These 
natural parameter values thus give an appealing expla- 
nation of both the smallness of neutrino masses and the 

Figure 7: 

suppression of LFV decays. G 2 (blue) 



masses 



Mi 



to , 




The functions ^(red) and 



4.2 LFV Decays k -> /// 

Another class of processes that are of interest for our model are rare flavor violating decays 
of the type \x — > eee. As in the case of the processes /$ — > fjj the allowed decay channels 
are restricted by the flavor symmetry. If we do not consider the heavily suppressed diagrams 
that couple to VEVs in the neutrino sector, it is clear that the process \i — > eee is not allowed 
by the Z3 symmetry of the charged lepton sector and the most constraining process is given 
by — >■ fi~e + . 

This process can be mediated at tree-level by the neutral components of <p" as depicted 
in Fig. 8a and its branching ratio is given by [22; 24] 



Br(r" 



3Qm 2 m 2 \ 

" 1 Brfr 



Ml 



fiw) = 1.7. 10~ 8 



62 GeV 



Mi 







(4.12) 



where we have used Br(r —> \ivv) 
1.7- 10" 8 [341, the effective mass 



0.174. Compared to the experimental upper bound of 



sin 2 a 



mi 



7 In 



A x a = was assumed, which implies a — n/4. 



cos 2 a 



TO 



$ 2 



(4.13) 
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(a) Tree level contribution of cp" (b) One loop contribution of S-r) 

Figure 8: Lepton Flavor Violating rare decay t~ — > fj,~ /i~e + . 



is thus only weakly constrained. All other processes mediated by ip'^ are further suppressed by 
y e y T or y^y & . Rare LFV processes mediated by these fields are therefore naturally suppressed 
by smallish Yukawa couplings and do not put a serious constraint on the model. 
Let us also estimate the magnitude of the diagram in Fig. 8b mediating r — > fJ-^ie, as this 
diagram may in principle be larger because it is not suppressed by Yukawa couplings that are 
known to be small. 

To get an estimate, we work in the limit of degenerate r\ masses M\ = M 2 = M3 = m,o 
and find 

9 3 

K Z. U* U I 



1 h h * h h * H(M s /m ) 

L n rijn tiik n e i k n i 

j t k=i i,i=i 



16vr 2 



where H(Ms/mo) is a dimensionless loop integral and 



Br(r — > fi \i e + ) = Br(r — > \wv 



T(t- 



rr 



Evaluating the sum, we find X^,fc=i YLi,i=\ ^rijh^h^h*^ = {h\ - h\h\ + hf) and the 
experimental bound 



Br(r" 



H \i e 



159 GeV V(^ + hi-h 2 1 h 2 2 )H(M s /m ) 



m 



1.7-10" 



(4.14) 



can easily be evaded even for small values of mo ~ 178 GeV ~ 308/\/3 GeV (which would give 
the correct dark matter relic abundance of r/ in the degenerate limit we are considering here, 
as will be discussed in sec. 5.2) and order one Yukawas (assuming H(Ms/mo) ~ 1). For the 
parameter ranges preferred by one-loop neutrino mass generation, i.e. hi ~ 0.1, the expected 
branching ratio is too small to expect a signal in next-generation experiments. In summary, 
we can conclude that the flavor symmetry effectively protects against lepton flavor violating 
interactions. 



4.3 Anomalous Magnetic Moment of Muon 

Let us now briefly discuss the anomalous magnetic moment of the muon. The contribution 
from the exchange of the neutral component of ip" should give the largest contributions, as 
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it is proportional to the tau Yukawa coupling squared. It has been calculated previously [IT] 
and amounts to 



Aa, 



2V2- 



7T 



Ml 



1.5 • 10 



-12 



100 GeV\ : 



Ma 



(4.15) 



which is negligible and cannot account for the reported deviation of (290 ±90) x 10 -11 [46; 47], 
from the Standard Model. The charged components of rj also contribute to the anomalous 
magnetic moment of the muon, with a strength given by [43; 47] 



m: 



3(4vr 



2^2 



h 



Ml 



\mJ + M| 2 



1 F 2 



-1.8 x 10 -12 ^ 



hj 
0.1 



100 GeV\ 
Mi ) 



2 f F „(Ms- 
r2 ^Mi ■ 



F 2 (l) 



(4.16) 



This therefore gives a very mild constraint on the masses and Yukawa couplings of the 77's. 
In the preferred parameter space for neutrino mass generation, this contribution is negligible. 
Note that the contribution goes in the opposite direction of the reported excess and it can 
therefore not be used to explain it [43]. 



5 Dark Matter 

In this section we discuss dark matter candidates of the model and their phenomenology. 



5.1 Dark Matter Candidates and their Stability 

To start off the discussion of possible dark matter candidates in our model, let us dwell on 
the remnant symmetries left over after symmetry breakdown. While the Qs x part of the 
symmetry group is completely broken 8 , there is a Z2 symmetry given by 

TZ: L^-L t c ^-t r\i -»■ -T]i , (5.1) 

which is the (— 1) L ' remnant of the auxiliary Z4 symmetry i L , where V = L + N„ is the 
generalized lepton number symmetry that is the sum of the usual SM lepton number with 
the r\ number N„. At the renormalizable level after symmetry breaking, there is another Z2 
symmetry of the model given by 

A : S -> -S rn-> -rji ■ (5.2) 

This is purely an accidental symmetry that emerges due to the particle content and the 
requirement of renormalizability and not a remnant of some symmetry we have imposed on 
the model. The reason why it emerges can be traced back to the fact that the SM fermions as 
well as x transform only under the generators S and T that form the subgroup and thus 

8 There have been several studies of dark matter, which is stabilized by a remnant subgroup of a flavor 
symmetry [48], while it is completely broken in our model. 



19 



there are no operators of the form f>0^ , where ip is a field transforming non-trivially under 
X (e.g. fields transforming as 3j with i ^ 1 such as S and rj) and Oa± is an arbitrary operator 
formed by fields transforming under A±. The symmetry A makes the lightest component of 
S and rj stable, which implies that the dark matter candidate is either fermionic or bosonic. 
This symmetry, however, is only an accidental symmetry and there is thus no reason for higher 
dimensional operators to respect this symmetry. Such a higher dimensional operator O with 
A[0] 7^ O would lead to a decay of the dark matter candidate. On the contrary, all higher 
dimensional operators have to respect the symmetry 1Z[(D] = O, as this symmetry is a remnant 
of an exact symmetry and is therefore also exact. We will now show that this requirement 
pushes up the dimensionality of the higher dimensional decay operators to a level where the 
dark matter candidate is stable for all practical purposes. Since the discussion depends on 
whether the dark matter candidate stems from r\ or from S, we discuss the two possibilities 
in turn. 

Scalar DM: Any effective operator that would mediate a decay of the lightest component 
of rji has to be of the form 

= m O§kf l {0 4 , k 4> l ) (5-3) 

where (0 ( p k( j, l ) is built out of SM-singlet flavon fields and transforms even under 1Z. As rj is 
odd under TZ, the operator 0^ =1 , which is built up of SM fields, has to be also odd under 
1Z to make the complete operator invariant. Obviously the complete operator O is odd under 
the accidental symmetry A and thus mediates DM decay. 

Since 1Z acts upon SM particles as the discrete subgroup of lepton number (— 1) L , the 
operator 0^^ =1 has to violate lepton number by an odd unit and has to transform as an 
electroweak doublet. The lowest dimensional operators in the SM arise at dimension six and 
violate L by one unit (See [49] for a recent review of gauge invariant dimension 6 operators.) 

Lu c d c d c Ld c d c d c LQQd c e c Qd c d c (5.4a) 
X ] LQQQ X ] e c u c u c d c X ] LQu c d c x ] e c QQu c X ] LQu c d c . (5.4b) 

All dimension 6 operators in Eq. (5.4a) break baryon number by one unit, B — L by two units 
and preserve B + L. The dimension 7 operators in Eq. (5.4b) on the other hand break baryon 
number by one unit, preserve B — L and break B + L by two units. They are formed by 
adjoining x to a dimension 6 proton decay operator. Since baryon number is an accidental 
symmetry in our model (in the same way as in the Standard Model), these operators are never 
generated 9 within the model and thus dark matter is stable within the model. They rather 
parameterize some baryon number violating physics, which from proton decay experiments is 
pushed to scales of the order of As « 10 16 GeV. 

To form a singlet under the flavor symmetry, the second operator is needed to make 

the total operator O a singlet under the flavor symmetry, as rji transforms under X while 
C^ =1 does not. It has to be composed of an even number of flavons as under the Z2 
subgroup generated by 10 X 2 only transforms non-trivially. 

9 Except through instantons and sphalerons, which do not play a role here, in the same way as in the SM. 
10 This element generates the center of the group and thus commutes with all group elements. 
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If we assume the presence of baryon number violating operators at scale A#, the dark 
matter candidate r\ decays into quarks and one lepton. Under the assumption that the flavor 
part of the operator is related to the breaking of the flavor symmetry Ap, a DM decay operator 
formed by a dimension 6 SM operator O^f 1 is suppressed by A^: 



Hence, the lifetime of DM can be estimated to be 



(5.5) 



(M' = 1.9. 10- GyI ( *lA* ] (5.6) 



m 7 v \{<i> k <t>i)) J \W 16 GeVj \ m v J \(c(> k< j,, 

and the dark matter candidate is thus stable even on cosmological time-scales, if one assumes 
'traditional' values for the scale of baryon number violating physics. However the operators 
in Eq. (5.4a) are not those directly tested in proton decay experiments and the physics of 
baryon number violation might be such that the operators in Eq. (5.4a) are suppressed by 
a smaller energy scale than the one responsible for baryon decay. We will come back to the 
issue of induced proton decay at the end of the subsection, but now we want to turn the logic 
around and derive bounds on Ag and Ap from the fact that dark matter is still around. 

Decaying DM models are constrained by WMAP to r 1 > 123 Gyr at 68% C.L. [50] 
and WMAP+SN la to > 700 Gyr at 95.5% C.L. [51]. Furthermore, decaying DM is 
constrained by possible neutrino final states [52], which serve as a conservative limit, since 
neutrinos are the least detectable SM particles. The exact bound depends on the DM mass 
ranging from 10 22 s = 10 8 Gyr at 0(1 GeV) and increasing almost linearly on a log-log plot to 
10 28 s ~ 10 14 Gyr at 0(100 TeV). Diffuse gamma ray constraints from Fermi data yield a limit 
of r _1 > 10 26 s ~ 10 12 Gyr [53] for the decay into a pair of charged leptons. Here, DM decays 
into one lepton and quarks, which might lead to further softer leptons in the final state. Hence 
the bounds do not directly apply, but we will use it to obtain an order of magnitude estimate 
for the suppression scale of the lowest order DM decay operator in Eq. (5.3). Using the limit 
from diffuse gamma rays with r^ 1 > 10 26 s benchmark value, we obtain a limit on the 
suppression scale of 

1/5 



[A B A F ) >6 10 GeV^p^j l v( i 00 GeV) 2 J 



(5.7) 



Because of the high dimensionality of the operator, the bound on the suppression scale Ab,f 
does not depend strongly on the bound on the lifetime. 

All of the operators in Eq. (5.4) lead to DM induced proton decay 11 into a final state 
lepton and final state mesons 

r]i + N^L + M. (5.8) 

As the proton as well as the DM are non-relativistic and they annihilate at rest, the induced 
proton decay leads to similar kinematics as in the ordinary proton decay, but the total rest 



11 Induced proton decay has been studied in the context of asymmetric DM [54]. However, their analysis 
does not apply in our case, because the induced proton decay is mediated via a different operator with different 
kinematics, since one of the final state particles has a non-negligible mass of the order of the proton mass. 



21 



energy E ~ m„ + tun ~ rn v is much larger compared to the ordinary proton decay with 
E ~ m^v- Hence, the final state particles appear to originate from the decay of a much 
heavier particle and the experimental signatures change. Therefore, the existing limits on 
proton decay are not directly applicable. However, in generic GUT models, for example, the 
operators given in Eqs. (5.4a), (5.4b) and the proton decay operators are generated at the 
same energy scale. 

Fermionic DM: Similarly to scalar DM consisting of the lightest component of rji, S can 
decay via higher-dimensional operators. They are generally of the form 

SOsm (O^) , (5.9) 

where Osm transforms like a spin | fermion, which is a singlet under the SM group, but 
transforms non-trivially under the flavor symmetry 12 . The lowest dimensional operators Osm 
emerge at dimension | 

u c d c d c QQd c xQu c d c xQQQ- (5.10) 

Note that these operators transform trivially under 1Z, as does S. All of these operators violate 
baryon number by one unit and therefore, they lead to induced proton decay. However, the 
kinematics is quite different compared to ordinary proton decay, because the lowest order 
operators do not contain a final state lepton. 

Similarly to the scalar case, there are bounds from astrophysical observations. As DM 
decay only arises at dimension 8, the bound on the suppression scale does not depend strongly 
on the exact bound on the lifetime. Therefore, we again make a rough estimate of the bound 
on the suppression scale by using the same lifetime as in the scalar case and we obtain 

(A^)-> 9 ao 8G ev(^)-( I JM^) 1/1 (,n, 

due to the lower dimensionality of the DM decay operator. 
5.2 Dark Matter Phenomenology 

We now give a brief overview of the phenomenology of the two different dark matter can- 
didates. We will estimate the DM abundance and detection possibilities for the different 
scenarios and show that there is a region of parameter space where the correct abundance 
can be obtained. A detailed calculation is beyond the scope of the present work. Again, we 
discuss the different dark matter candidates separately. 

Scalar DM: 

The scalar dark matter candidate is a component of an inert EW doublet. Therefore, we 
are going to translate the analysis for scalar multiplet DM done in [55] to our setup. A detailed 

12 Note that S transforms under the symmetry generator X, while Osm does not. Therefore the operator 
(CVfc^,) is needed to form a singlet. 
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analysis would require the precise calculation of the rji mass matrices. We assume that one 
of the triplets rji is sufficiently lighter than the other two, such that we do not have to take 
them into account during freeze-out of DM, i.e. they have to be at least 20% heavier than the 
DM candidate [56]. In the following, we will denote the triplet containing the DM candidate 
by r)DM with direct mass term M VDM . We are going to assume, as we did previously in the 
section about the neutrino masses, that the direct mass term M r)DM dominates over all mass 
terms induced by VEVs. Hence, the DM mass is approximately given by the direct mass term 
M r]DM . In the limit that the mass splittings are below 1%, we can neglect the annihilations 
via other scalars and concentrate on the pure gauge (co) annihilation channels. Following [55], 
there is an upper bound on the DM mass of an inert doublet of m* = 534 ± 25 GeV(3cr) from 
overdosing the universe in this limit. The correct DM abundance is obtained for m* . As t]dm 
is in a triplet representation of Qs x A4, there are three almost degenerate doublets, which 
all contribute to the DM density equally. Therefore, the upper bound on the DM mass is 
lowered by approximately a factor of \/3 to m* ~ 308 GeV, which is consistent with direct 
searches for scalar particles, as discussed in sec. 7. 

Today, the mass splitting between DM and the next-to lightest particles forbids gauge 
interactions kinematically due to the small DM velocities, unless it is tuned to be very small 
(< 0(100) keV), and DM can only be detected via the couplings to scalars, specifically via 
the Higgs portal. The spin-independent cross section for scattering of DM off the neutron is 
given by [57] 

(5.12) 

with Xl being the coupling of DM to the Higgs, \x the reduced mass of the DM- neutron system, 
m p the mass of the nucleon, ran the mass of the Higgs and / parametrises the nuclear matrix 
element, 0.14 </ < 0.66, which we took from [57]. The estimated cross section is well below 
the current experimental limits by XENON100 [58], which is the most sensitive DM direct 
detection experiment in this mass region. 

Note, the discussed parameter point is only an example that proves the possibility of 
obtaining the correct DM relic density. For larger mass splittings, the annihilation via scalar 
interactions cannot be neglected in the calculation of the DM relic abundance and the direct 
detection cross section is enhanced. 

Fermionic DM: For the discussion of the fermionic DM candidate contained in S, we 
follow the discussion in [19] to show that it is possible to obtain the correct relic abundance. 
For completeness, we repeat the relevant steps with the necessary changes. At tree-level, 
there is only the mass term \/3MsSS = Ms(S± + + S%) and thus all components of 
S are degenerate. At loop-level this degeneracy is lifted and for concreteness we here take 
> M<j 2 > M§ , where Si are mass eigenstates. The states $2,3 can decay into Si and 
leptons by the interchange of r\ and thus at the present time only Si is around. However, 
due to the near degeneracy, the freeze-out of all three species runs in parallel. Coannihilation 
processes of the type SiSj — > SM with i 7^ j are suppressed in comparison to annihiliation 
processes SiSi — > SM, because they require an additional mass insertion along the 77 line. It 
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is thus a very good approximation to consider the freeze-out of each component separately 
and the total relic abundance is thus just given by the sum of the abundances of S±, S2 and 
S3. 

The annihilation cross section for each Sk into leptons in the limit of vanishing lepton 
masses and scalar mass splittings [59] is given by 

;i-2r; + 2rf) 



(av) = bv z + 0(v 4 ), 
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In the limit of Ms <C M%, the expression for the p-wave simplifies to 
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i.e. the cross section scales with (hi/rrii) 4 . The relic density of the SM singlets S, taking into 
account the mass degeneracy of the components of S, can then be obtained from [60] 

n° s M s . 2 



with n s being the number density of S today, which is 



r 1 = ( 



(3< 



0.088 gl /2 M Pl M s 3b 
xjs 



(5.15) 



(5.16) 



where so = 2970/cm 3 is today's entropy density, the critical density is p c = 3H 2 /(8ttG) = 
1.05 • 10~ 5 /i 2 GeV/cm 3 , the Planck mass M P[ = 1.22 • 10 19 GeV and the dimensionless Hubble 
parameter h. At the freeze-out temperature, the ratio Xf = M$/T is determined by 

0.0764M«(66/^/)c(2 + c)M s 
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with the effective number of degrees of freedom g* at freeze-out. After eliminating the cross 
section with Eq. (5.16) and Eq. (5.15), we obtain 
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Following the discussion in [19; 60], we rewrite Eq. (5.15) and Eq. (5.18) as 
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10 and c = 1/2. We solve these equations numerically for fixed values of h\ = J12 
and Mi = M2 and show the resulting contour lines with the correct DM relic abundance 
in the plane M\jh\ = M2//12 vs. Ms in Fig. 9. Hence, it is possible to obtain the correct 
DM relic abundance for fermionic DM, although large Yukawa couplings hi are required. 
Similarly to the scalar DM scenario, we expect the cross section to raise with non-vanishing 
mass splittings of the scalars r/i, which allows for smaller Yukawa couplings hi. 
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6 Extension to Quark Sector 

So far we restricted ourselves to the discussion of the flavor structure in the lepton sector. 
Given the different structures in the lepton and quark sector, one might wonder whether 
and how this model can be extended to the quark sector. In the following, we will discuss 
a few simple possibilities to incorporate the quark sector without enlarging the flavor group. 
It is necessary to specify how quarks transform under the flavor symmetry as this will to a 
certain extent determine the collider signatures of the model. Alternatively, it is interesting 
to look for a group extension of the flavor group, which preserves the structure in the lepton 
sector, but allows for new structure in the quark sector [32]. Here, a viable extension of the 
full flavor group Q 8 x A 4 is the group Q 8 x T = SG(192, 1022) [32] being the analogue of 
the extension of to T', which has been used to explain the flavor structure of quarks and 
leptons simultaneously [61]. A detailed discussion of quark flavor observables is postponed to 
future work. 

Quark Sector Mirroring the Lepton Sector 

We can use the same assignment for the quarks as for the leptons with respect to (Qs x A4) x Z4, 
i.e. 

Q~(3 1j 1), u c + c c + t c ~ (1 1 + 1 2 + 1 3 ,1) , d c + s c + b c ~ + l 2 + l3,l) • (6.1) 
This assignment leads to the following Yukawa couplings in the Lagrangian 

- C q = y u Qxu c + y c Qxc c + Vt Qxt c + Vd Qxd c + y s Qxs c + y b Qxb c + h.c. , (6.2) 
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which amount to the mass matrices of the quarks 



V2 



fiydiag(y u , y c , y t ) 



and 
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(6.3) 



Hence there is no mixing in the quark sector, i.e. the CKM mixing matrix Vckm = Vj V u = 1, 
which is a good leading order approximation to the CKM mixing. The Cabibbo angle can 
be produced by a cross-talk of operators from the neutrino sector [29], e.g. the operator 



(Qxd c 



6 2 ) 2 /M 4 leads to a non-vanishing Cabibbo mixing angle. It has to be of the 



order ((frifo) 2 / (M 4 ) ~ 10~ 4 (m s /95 MeV), in order to generate a large enough mixing in the 
down-type quark sector to explain the Cabibbo angle. Within the model, the operator can 
be generated at one loop with <//W running in the loop. However, the contribution turns out 
to be too small and a different mechanism is required to generate this operator. 

Flavor changing neutral currents are naturally suppressed at the leading order, since there 
is a selection rule ADASAB = ±2 as well as AUACAT = ±2 in the flavor basis for four 
Fermi operators similarly to the lepton sector. It has been claimed in [22] that leptonic Kaon 
decays result in a relatively strong bound of Mq > 510 GeV on the effective mass Mq defined 
in Eq. (4.13). However, there is an error in the calculation. The final result should not depend 
on the Kaon mass mx but and the corrected expression in our model reads 
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The branching fraction is constrained to be less than 4.7 -10 . Using m s = 95 MeV, 
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Quarks Transforming under generator X 

Another interesting possibility that is not possible in A4 models is to assign the quarks to 
representations that also transform under the group generator X. Since the top mass is large, 
we want it to be generated at the renormalizable level, while all the other quark masses might 
well be the result of higher order effects. Looking at the multiplication rule 



3j x 3j 



fc=i 



(6.6) 



it is clear that if one assigns Q ~ (3_2, 1) and U c ~ (33, 1) there is only one Yukawa coupling 
at the renormalizable level 

-C t = y t QxU c + h.c. , (6.7) 

which generates the top mass. The charm and up quark masses, as well as up sector mixing 
are generated by operators of the form 

-C u = [Q X U c {Ml)\i + vt ,2) [QxU c {h<t>2)]i + h.c. , (6.8) 
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where the sum goes over all singlet contractions of the fields. There are certainly enough para- 
meters to fit the quark masses and up-type mixing. Actually, there are no further predictions 
besides the large top mass, since there are too many free parameters. 

In the down-type sector we can either utilize the same structure as in the up-type sector 
or, as the bottom quark mass is closer to the charm mass than to the top mass, we can use 
the assignment D c ~ (3]_, 1). With this choice there is no tree-level operator of type (6.7) 
allowed and all down type quark masses and mixing arise from 

-C u = yf^ [QxD^Ml)^ + vf ' 2) [QxD c (hh)]i + h.c. (6.9) 

We will not discuss this possibility further here, as we are primarily focused on the lepton 
sector. 

Additional EW Higgs Doublet H q ~ l x 

Another possibility is that the flavor structure in the quark sector could be completely unre- 
lated to the one in the lepton sector. In particular, the quarks might not transform under the 
flavor symmetry in the lepton sector. This can be simply achieved by assigning the quarks to 
the singlet representation of the flavor group. In order to generate the quark mass matrices, 
we have to introduce an additional EW Higgs Doublet H q , which does not transform under 
the flavor group. Hence, the flavor structure in the quark sector is unchanged compared to 
the SM one. Therefore, we do not discuss this possibility further and we will only briefly 
comment on its collider phenomenology in sec. 7. 

The only effect 13 of the additional Higgs doublet on the discussion in the preceding sections 
is to rescale the VEV of H such that 

(H°) 2 + (H°) 2 = \{^2G F )^ = ^(246 GeV) 2 

is maintained. 

7 Collider Phenomenology 

Our model predicts several new particles with EW charges at the EW scale. In this section, 
we will concentrate on the simplest extension to the quark sector given in sec. 6, where quark 
doublets are assigned to the triplet representation 3^ of the flavor group and obtain their 
masses from a coupling to the flavored Higgs x, as discussed in the previous section. We will 
briefly comment on the possibility to have a separate Higgs for the quark sector in the sec. 7.5. 
Besides the fermionic singlets S, there are several EW doublets, which can be grouped in three 
different categories, the Higgs h, which obtains a VEV, the two partners of the Higgs in the 
flavor triplet x, namely <p' and (p" , and the additional inert EW scalar doublets r/. In the 
following, we sketch the different production and decay channels and discuss their implications 

13 Here we assume that H q does not give a leading order contribution to the Weinberg operator. Symmetries 
can always be adjusted in order for this to be the case. If H q does give such a contribution there will be one 
more free physical phase in the neutrino mass matrix that cannot be rotated away. 
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for direct searches at colliders as well as the current bounds on the existence of new particles 
beyond the SM. However, a detailed study is beyond the scope of this presentation. 

After a brief discussion of electroweak precision constraints and a short summary of the 
main experimental results, we will discuss each class of new particles separately. 

7.1 Electroweak Precision Constraints 

The experimentally measured values of the oblique parameters S and T have been obtained 
by several precision measurements at LEP and Tevatron. The PDG [34] quotes values of 
S = -0.04 ± 0.09 and T = 0.07 ± 0.08 at 95% C.L. with a correlation between S and T of 
88% for a reference value of rrih,ref = 117 GeV. 

A general discussion in a multi-Higgs doublet model with an arbitrary number of Higgs 
doublets with hypercharge Y = ±i and an arbitrary number of SM singlets has been given in 
[62]. The expressions for the oblique parameters can be directly applied to this model, since 
the flavor symmetry only leads to additional restrictions on the masses and mixing matrices. 
We only estimate the contribution to S and T in the limit of small mixing in of f], fy, H and 
only consider the mixing of if' with (p", which exactly corresponds to the region in parameter 
space being studied in the previous sections. In this limit also the charged and neutral scalar 
masses of the doublets f] coincide. In this approximation, the contribution of H exactly 
cancels with the subtracted SM term, the contribution of fj and (pi to T vanishes, and fa 
does not contribute to S, since it does not couple to the gauge bosons in this approximation. 
Hence, the final contribution to S originates from fj and is given by 



where a = 1, . . . n sums over the EW doublets contained in fj with the charged scalar masses 




(7.1) 



m a and 9w denotes the Weinberg angle. The function G is defined by 




(7.2) 



Its absolute value is monotonously decreasing for z — > starting from (5(1) = —0.216 to 
(5(0) = 0. The contribution from ip'W to S and T is given by 




(7.3) 




(7.4) 



where the mixing matrix in the neutral y/W sector, is defined in Eq. (2.7) and m a (fib) 
denotes the charged (neutral) masses of the fields contained in y/W . The function F is defined 
by 




(7.5) 
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The next-to leading order corrections are suppressed by small mixing angles in the scalar sec- 
tor. Hence, the model is consistent with electroweak precision tests in the phenomenologically 
interesting region, i.e. for small mixing in the scalar sector. 

7.2 Summary of Relevant Experimental Results from Colliders 

Recently, after the initial announcement of a SM-Higgs like resonance by ATLAS [ I] and 
CMS [ ], which was mainly based on the diphoton as well as the h — > ZZ* — > 41 channel, 
several other channels have been measured or updated [4; 6]. We will briefly summarize the 
current status. The current best fit values for the mass of the resonance are 126.0±0.4(stat.)± 

0. 4(sys.) GeV by ATLAS [3; 4] and 125.8 ± 0.4(stat.) ± 0.4(sys.) GeV by CMS [6]. The results 
are usually reported in terms of the signal strength normalized to the SM prediction, i.e. 

_ a(pp — > h)Br(h — > X) 
a(pp -> h S u)Bi(hsM -> A) ' 

The two main channels are the decay into two photons and h — > ZZ* — > 41. While the 
h -» ZZ* -)■ 41 rate seems to agree with the SM prediction with R zz = 1.2 ±0.6 for ATLAS [ ] 
and Rzz = 0.8^Q2g for CMS [6], the h — )■ 77 rate seems to be enhanced with i? 77 = 1.8 ±0.5 
for ATLAS and i? 77 = 1.56 ± 0.43 for CMS. The remaining channels include h -)• WW* with 
a signal strength of R ww = 1.4 ± 0.6 (ATLAS) and R ww = 0.74 ± 0.25 (CMS) and the two 
channels with decays into fermions hV — > bbV with a signal strength of R^y = —0.4 ±1.1 
(ATLAS) and R bW = 1.3±g;| (CMS) as well as h -> rr with R Tfjj = 0.7 ± 0.7 (ATLAS) and 
R T fjj = 0.72 ±0.52 (CMS). All channels but the decay of the Higgs boson to two photons are 
in agreement with the SM prediction. The deviation in the diphoton channel is intriguing, 
as in the SM this decay proceeds via a loop diagram and is thus sensitive to new physics 
contributions. However, so far, the deviation is at the 1 — 2a level [63-65], if the uncertainties 
are taken into account conservatively. Besides the discovery of a Higgs- like resonance, the 
LHC has put strong constraints on any physics beyond the SM. 

Charged Higgs particles are constrained by searches at LEP and LHC. At LEP, charged 
Higgs particles are produced via a virtual Z* in the s-channel, i.e. e + e~ — > Z* — > H + H~ , 
and studied via their decays into tv t as well as cs assuming their branching ratios add up to 

1, i.e. Br( H + t + v t )+Br( H + -> cs )=1. This results in a bound of m H + > 79.3 GeV [34]. 
Independent of any assumptions on the branching ratio, the invisible Z decay leads to > 
45 GeV [34]. CMS searched for charged Higgs particles [66], which are produced in top decays, 
t —> H + b and constrains their branching ratio Br(t — > H + b) to less than 2%-4% for charged 
Higgs masses between 80 and 160 GeV. Similarly, the search by the ATLAS experiment [67] 
yields bounds on the branching ratio Br(t — > H + b) of the order of l%-5% for charged Higgs 
masses in the range between 90 and 160 GeV, assuming Br(i? + — > t + v t )=1. 

7.3 EW Higgs Doublet H 

We will first consider the limit in which there is no mixing between the Higgs h and the flavons 
(pi. In the limit of no mixing, the tree-level couplings of the Higgs h contained in the EW 
Higgs doublet H to gauge bosons are identical to the SM couplings. In addition, the flavor 
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conserving tree level couplings of the Higgs h to fermions also agree with the SM ones. Note 
that there might be small corrections, since quark mixing vanishes at leading order and the 
Higgs couplings conserve all flavor numbers separately. As there are no new colored particles 
and the coupling of the Higgs to ti is the same as in the SM, the loop-induced coupling of the 
Higgs h to gluons agrees with the SM one. In summary, the production of the Higgs h as well 
as all tree- level decay channels and the decay into gluons are exactly like those in the SM. The 
Higgs decay into two photons is the only decay channel that can show a significant deviation 
from the SM in this approximation. If any of the other new scalars were light enough, there 
would be additional tree level Higgs decays into pairs of these scalars and such scenarios are 
therefore constrained. The decay h — > SS , if kinematically allowed, is loop suppressed. 

Mixing of the Higgs h with the flavons (pi leads to a suppression of all tree-level couplings 
to gauge bosons and fermions. Hence, the production cross section is reduced according to 
the admixture of the flavons 4>i to the Higgs boson. As Higgs decays into ZZ* are close to 
the SM value, the admixture of the flavons <pi to the Higgs h is limited. 

Finally, let us discuss the diphoton decay channel. The SM contribution is dominated by 
the W boson contribution and the smaller top loop contribution, which interfere destructively. 
In our model, the decay into two photons receives additional contributions from charged 
scalars in the loop, which are contained in the EW doublets cp', ip" as well as 17. Any enhancing 
contribution has to interfere constructively with the SM W boson loop or dominate over the 
W boson contribution. The contribution of additional charged scalars pi with a charge one, 
coupled to the Higgs boson via the Higgs portal 

Pi = c p ^H\p t \ 2 , (7.6) 

has recently been studied in [68]. The ratio of the effective coupling of the Higgs boson to 
two photons vs. the SM prediction is given by 



1 -/Z c pM m P t ', 



(7.7) 



where the function h is depicted in Fig. 10. To obtain an enhancement of a factor of 2 (1.5), 
one thus needs a value of 

E, , . —0.41 (—0.22) for constructive interference I ,„ „. 

c p .h(m p ) = < . , > . (7.8) 

. I 2.41 (2.22) for destructive interference j 

Hence, a large negative coupling c p ~ —2 is necessary to obtain an enhancement factor of 2 
for a single singly charged scalar of mass 100 GeV. Such a large negative coupling destabilizes 
the vacuum and leads to charge breaking minima unless \c p \ < \fW p ~ ^J~\ p j1 is fulfilled, 
where X p denotes the quartic coupling A p |/o| 4 /2. Note that this requires very large values for 
A,. 

Let us now use this formula to estimate the deviations from i? 77 = 1 that can be expected 
in this model. In total we have 11 charged scalars, 9 from the doublets f\ and two from the 
doublets ip' , ip" . The interaction of the last two scalars with the Higgs field can be expressed 
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Figure 10: Plot of function h of Eq. (7.8). 
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Figure 11: i? 77 in the case where all charged scalars 
have the same common mass M as a function of 
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with m ,(/) + defined in Eq. (2.6). In the limit of large m 2 ,^ these two fields contribute 



cy/i(Mi) + c^f'hiMt) = 0.1 + 



22GeV\ 2 /22GeV\ 2 
M% J + \ M°l J ' 



(7.10) 



The couplings of the charged components of the rj fields are given by 

3 

W|2 



(7.11) 



i,J=l 



as dictated by the symmetry. The coefficients are essentially unconstrained except for 
the fact that the combination that couples to the DM particle should not be too large, to 
avoid the bound from direct detection. In the limit where all charged scalars have a common 
mass M, we see from Fig. 11 that M = 200 GeV requires | Yl j=i c vj = — 1-46 (—0.95) for 



R 



77 



2(1.5). 



In case the h — > 77 anomaly persists, it would be interesting to measure h —> ~fZ, since 
it originates from similar diagrams, where one photon is replaced by one Z boson. A cross- 
correlation of the two measurements would allow one to determine the isospin of these par- 
ticles. In our model all charged scalars are part of SU(2) doublets allowing us to distinguish 
it from other models which have EW multiplets in the loop with different EW charges, like 
singlets or triplets. 
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7.4 Further Scalars 



Besides the Higgs h, there are several additional scalars, such as the flavor-violating EW 
scalar doublets ip'^ as well as fj, which do not acquire a VEV, and the flavons fa, which 
acquire a VEV. See appendix A for the scalar mass spectrum. 

Flavor- violating Higgs Doublets ip r ( r > : The neutral components of the flavor- violating 
Higgs doublets ip'V' have neither tree-level couplings to tt nor couple to two EW gauge bosons 
at tree level. Hence, they are not produced by any of the standard Higgs production channels, 
but they can be produced via associate production with two different quarks, gg p'^qq' . 
The dominant channel is gg — > (p'iu, <p"ic, which has a cross section of the same magnitude as 
production of a Higgs boson in association with a tt pair. Therefore, they are not constrained 
by the current heavy Higgs searches. Other productions channels are qq' — > p^'\ where q and 
q' are in different generations as well as pair-production in vector boson fusion WW, ZZ —> 
ip'(p',(p"ip". These processes are suppressed compared to the main Higgs production channels 
at the LHC. Note, however, that the decays of the flavor-violating Higgs doublets y/W might 
lead to distinct flavor-violating signatures similar to the recent analyses of flavor-violating 
Higgs decays in models with flavor symmetries [25-28; 69]. 

As the charged Higgs particles contained in tp'w do not couple to tb, the LHC limits do 
not apply. Hence, the charged Higgs particles in our model are only constrained by the LEP 
limits discussed previously. 

Although there are no constraints yet, upcoming searches will test the allowed range of 
masses, because the flavor-violating Higgs doublets (p'^ stem from the same flavor triplet 
as the Higgs doublet H, and therefore their masses are determined to be given by scalar 
couplings times the EW VEV. Their masses may therefore not be raised arbitrarily high, as 
discussed below Eq. (2.6) 14 . 

EW Scalar Doublets r): The neutral components of the EW scalar doublets fj do not 
couple to quarks and particularly not to tt as well as two EW gauge bosons. They can be 
pair- produced in vector boson fusion WW, ZZ — > r)ir\i. Hence, similarly to the flavor- violating 
Higgs doublets, they are not produced via the main Higgs production channels and the current 
bounds from heavy Higgs searches do not constrain fj. Also, the charged components of fj are 
not constrained by the current LHC searches, because they do not couple to quarks directly 
and the charged Higgs bounds do not apply. Therefore, they are only constrained by the LEP 
searches. 

Flavons 0^: The flavons fa do not have gauge interactions and they do not couple to 
fermions directly. However, they mix with the Higgs h, which is constrained by the Higgs 
searches to be small, since a large mixing suppresses the production cross section of h and 
therefore all rates relative to the SM expectation. In conclusion, the scalar mass eigenstates 
which are dominantly composed of the flavons fa are only produced via mixing with the Higgs 
h and thus there are no limits from current searches due small mixing. 

14 Note that if one introduces soft-breaking terms that respect the Z3 symmetry, it is possible to adjust the 
mass terms arbitrarily [22]. Alternatively one may introduce an EW singlet scalar that transforms as 3j and 
breaks to the same subgroup as X- This can be realized without introducing a vacuum alignment problem. 
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7.5 Variant with Additional EW Higgs Doublet H q 



As we discussed in sec. 6, another simple possibility to incorporate quarks in the model is 
by assigning all quarks to the trivial representation of the flavor group and introducing an 
additional EW Higgs doublet H q , which transforms trivially under the flavor group. This 
leads to different collider signatures compared to the previously discussed scenario. Soon, 
these scenarios can be experimentally distinguished at the LHC. We will highlight the most 
important differences. 

The discussion of the fermions S as well as the scalars fj remains the same. The main 
changes are in the Higgs phenomenology. In contrast to the other scenario, the component 
in x which obtains a VEV does not couple to quarks and therefore it is not produced in 
gluon fusion, unless there is mixing between x an d H q . Instead, the newly introduced Higgs 
H q will be produced in gluon fusion. In this setup, the observed resonance at 125 GeV 
would be associated with the mass eigenstate, which is dominantly composed of Higgs H q . 
As H q has exactly the same couplings to gauge bosons and quarks, but does not couple to 
leptons (especially r's), the decays into leptons are suppressed by the mixing between H q 
and H (contained in The diphoton branching ratio can be enhanced in the same way as 
discussed in sec. 7.3. 



7.6 Fermionic Singlets S 

The additional fermionic states S are SM singlets and only charged under the discrete flavor 
group. Furthermore, they only couple to lepton doublets and therefore their production 
cross section at hadron colliders is suppressed compared to colored particles and there are 
no relevant analyses at present. The production depends on the exact mass spectrum of fj 
as well as S. The production via t-channel fj exchange is always present in a lepton collider, 
e.g. e + e _ — > SS. If S is lighter than one of the components of fj, it is possible to produce 
S via EW production of these heavier components of fj and subsequently decay into S and 
one lepton. Unless S is the DM candidate, the fermionic singlet S will decay into a lepton 
and one of the lighter components of fj, which will subsequently cascade down to DM via EW 
gauge interactions. The signal is missing transverse energy and leptons (and possibly EW 
gauge bosons) in the final state. 

If S is lighter than all components of fj and therefore a DM candidate, there are bounds 
from mono-photon searches at LEP [ ]. As S only couples to leptons, the searches at 
hadron colliders are weaker due to the additional suppression from loops that couple leptons 
to quarks. The mono-photon searches at LEP probe the effective DM annihilation operator 
(eS)(eS)/A 2 , which are induced by the exchange of a scalar doublet rji^- The scale At of this 
operator is determined by A^ 2 = J^ fc \hk\ 2 /M% for Mj, 3> M$. The analysis in [ )] quotes a 
limit of (200 — 340) GeV for M$ < 90 GeV. Hence, this does not impose a strong constraint, 
since the smallness of neutrino masses points towards larger cutoff scales At . 



33 



8 Conclusions 



We presented a predictive renormalizable model of lepton flavor at the electroweak scale. 
The flavor group is extended in the scalar potential to Qs x A±, which allows a natural 
vacuum alignment at the EW scale [32]. This is the first model of its kind that explains the 
lepton flavor structure at the EW scale including the correct vacuum alignment. 

The SM Higgs boson is subsumed in a flavor triplet that couples to charged leptons (and 
quarks) at the renormalizable level, thereby eliminating the need to invoke higher dimensional 
operators, as is done in models with flavon singlets. Neutrino masses are generated at the 
one-loop level and are further suppressed by the fact that two small mass insertions are 
needed in the loop. This TeV seesaw is realized without imposing any new symmetries apart 
from the flavor symmetries. In the model there are five real free parameters, which gives a 
predictive framework and, in particular, a correlation between the atmospheric and reactor 
angle is predicted, which agrees well with the recent global fits. Furthermore the model 
automatically includes a WIMP dark matter candidate and its stability and phenomenology 
have been studied. It can explain the observed dark matter abundance and is consistent 
with current exclusion limits by dark matter detection experiments. Constraints from LFV 
experiments are loosened by the flavor symmetry in comparison to flavor generic multi-Higgs 
doublet models due to the remnant Z3 symmetry in the charged lepton sector. 

Finally, several possible extensions to the quark sector have been studied. We studied 
the collider phenomenology of the simplest extension to the quark sector, which does not 
require the introduction of new particles at leading order and in which the quarks multiplets 
transform like the lepton multiplets under the flavor symmetry, and commented on the other 
possibilities. We studied the possibility of the Higgs boson h to explain the observed resonance 
at 125 GeV, especially the enhanced diphoton rate, which can be straightforwardly explained 
by the multitude of additional charged particles contained in the EW scalar doublets, which 
all contribute to the radiative decay of h — > 77. The fact that the Higgs doublet is contained 
in a flavor triplet leads to distinct signatures at the LHC. There are additional EW scalar 
doublets (p'('\ which cannot be decoupled from the Higgs h, and therefore are accessible in 
searches at the LHC. As they do not acquire a VEV, they do not decay into gauge bosons, but 
only via Yukawa type interactions into fermions besides decays into other scalars. Because 
of the remnant Z3 flavor symmetry in the charged lepton sector, they only exhibit flavor 
violating decays into fermions in contrast to the Higgs h. 

It might be interesting to study leptogenesis in this model. Because of the flavor symmetry, 
the fermionic SM singlet S are degenerate in mass at tree level as well as all couplings but 
A2 are real. The degeneracy is only lifted at two loop order and therefore the induced mass 
splittings are small and there might be a resonant enhancement. This also introduces some 
CP violation into the mass matrix of S, but it has to be checked whether it is sufficient. We 
will leave a study of possible ways to obtain the baryon asymmetry of the Universe for future 
work. 
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A Vacuum Alignment and Scalar Spectrum of 
A.l Vacuum alignment 

The vacuum configuration given in Eq. (2.1) is naturally obtained from the most general 
potential 15 

V = P^(0i, 2 ) + V X ( X ) + v mbl (x, <k,<h) 
compatible with given symmetries, where Vx(x) is given in (2.4) and 

V^(01,02) =M?(010l)l + OLl (0101 )\ + ^ Oi(0l0l)3. " (01 01 ) 3- 

i=2,3 1 1 

+/X2(0 2 02)l 1 +/3l(0202)l 1 + V /?i(0202)3. "(0202)3. 

— 1 —1 1 1 

i=2,3 

+7i(0i0i)l 1 (0202)i 1 + ^2 7i(0i0i)3j " (0202)3. > 



i=2,3,4 



v x(x) =i4(xx)i 1 +pi(xxx)i 1 +Mxx)\ + Mxx)i 2 (xx)i 3 , 

V mix {x, 01, 02) = &3 (0101)1- {X ] X)l x + C23(0202)l 1 (X^l- ( AA ) 

compatible with given symmetries. The minimization conditions reduce to the equations 

a (o+ (a 2 + b 2 ) + a_ (a 2 - b 2 ) + 7+ (c 2 + d 2 ) + 7 _ (c 2 - d 2 ) + U x ) + Tbcd = 
b (a+ (a 2 + b 2 ) - a_ (a 2 - b 2 ) + 7+ (c 2 + d 2 ) - 7 _ (c 2 - d 2 ) + C/i) + Tacd = 
c (/3+ (c 2 + d 2 ) + /3_ (c 2 - (i 2 ) + 7+ (a 2 + 6 2 ) + 7 _ (a 2 - 6 2 ) + C/ 2 ) + Tabd = (A.2) 
d ((3+ (c 2 + d 2 ) - 13- (c 2 - d 2 ) + 7+ (a 2 + 6 2 ) - 7 _ (a 2 - 6 2 ) + U 2 ) + Tabc = 

v{M 2 + X x v 2 )=0 



with 



1 \/3 
Ui = ^ + ^Ci3V 2 for i = l,2 



15 We do not have to consider the part involving fj, because it does not change the minimization conditions 
of (j>i and Xi if it does not acquire a VEV. 
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and 



Mt 



2f4 + ( 13 (a 2 + b 2 ) + ( 2 3(c 2 + d 2 ), 



f - 6 P 



6 + 6 

2^3 ' 



7+ 



V37i + 74 
4^3 



IS 



72 + 73 j -pi 74 

ITT' and r= v!' 



with £ = a, f3. Since the number of equations matches the number of VEVs, vacuum align- 
ment is possible. Corrections to the scalar potential only arise on dimension 6 level. These 
corrections furthermore arise on one-loop level and are thus further suppressed. We therefore 
neglect VEV shifts arising from these interactions throughout this work. 

A. 2 Scalar Spectrum 

Scalar Spectrum — (f>i, x : Let us first discuss the visible sector, i.e. the flavons 4>i, fa, X that 
get VEVs and realize the symmetry breaking; the 77's are independent and will be discussed 
later. The fields can be classified according to remnant symmetries of the potential. There 
are the obvious symmetries 



^3 : X -> T 3 x, 



(A.3) 



with T3 = 0*rdiag(l, w 2 , uj)£$ t and 



Z 2 : 4>i -> S 4 4>i, x -> X, 



(A.4) 



with S4 = S7s 4 diag(l, 1, —1, — l)f2g 4 but there is another accidental symmetry of the potential 
V& not part of Qs x A4: 



Z 2 : 0i ->• O^i, 
with 16 O4 = S4 diag(l, 1, 1, -1)^ , where 



X ' X? 



It is useful to go to a basis 



9. 



5 4 



1 

71 



/ 1 -1 \ 
10 1 
-10 10 
\ 1 1 / 



(A.5) 



(A.6) 



(H, if',ip") T = & T x, (L e ,L fl ,L T ) T = ti! r L 



(A.7) 



where these symmetries are represented diagonally. Let us discuss the mass terms in turn: 



6 The alert reader will recognize this as an outer automorphism /14 denned in [71]. 
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the 9 physical scalars contained in x have been discussed following Eq. (2.6) Here we 
only report the expressions of the dimensionless couplings in terms of masses: 



2 



A xl2 = i ^3m 2 - 3^ m\ - 2m 2 m 2 + m\-A (M 2 - Ml) 2 + 3m 2 . - 2M 2 - 2M+J 
A x3l ,5 = -v / 3(M 2 +M%) 

A x3lA = -V3 (rn\ + \jm\- 2m 2 m 2 + m\-A {Ml - M 2 ) 2 + m\ - Ml - M%\ 
X xA = 6 (Ml - Ml) (Ai 

((pi)i and (02)4 transform as (1, —1, —1) and have a mass matrix given by 

/ 2(ac(v / 37A/-27 2 )+2fc72d) \ 

V . m l i((a,b,c,d,a 2 ) (c,d,a,b, (3 2 )) J 

with 

f- 2 ^ 2a 7A/ ( C -d)(c + rf) c r<A 6crd 

m u = -W3a a 2 + a^ (6 _ a)(ffl + fe) + -y j ~ ^ 

+ ^ (48v / 3a 2 6 2 - 3r (c 2 + d 2 ) + 8^72 (d 2 - c 2 )) 
(</>i)3 an d (</>2)3 transform as (1, —1, 1) and have a mass matrix given by 

/ 2(ac 72 -6d(72-2v / 37M)) \ 
" V3 ^ 

V . mn ((a, b, c, d, a 2 ) <-»• (c, d, a, 6, /3 2 )) / 



with 



2 26 2 (V3a 2 (a-6)(a + 6) + 2 7M (c-d)(c + d)) aaTd 

mn = 2V3a a 2 H ?. ; — 

b z — a z b 

2bcTd _ l_ r , 2 ^ + 72(c-d)(c + d) 

the real scalars /i, (4>i)i, (<Ai)2, (^2)1 and (^2)2 transform as (1, 1, 1) under the remnant 
symmetry. Here we don't give the full mass matrix but only give the mixing with the 
Higgs in the limit of small mixings. The mixing matrix with field / is given by 

tan 261/ = 2 J 2 (A.9) 



with 



_ _bvCl3 _ CwCl3 _ _dv(23 _ Cf ^23 

m h,(^i)i ~ ^ ' ^,(01)2 - ^ ' m fc,(<«i - ^/g ' m h,(Hh ~ ^ • 
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Scalar Spectrum — f): The relevant part of the scalar potential to calculate the mass 
insertions needed to calculate neutrino masses for the mass spectrum of t) has been given in 
Eqs. (3.5-3.8). To calculate the rj mass spectrum the complete interactions 

6V P = A i(x T c r 2'?x)l 1 (??r'720 ! %)i 1 + ^2e iax {x r cj2ax)^ 1 { r nla 2 Bm)\ 1 

+ AaC&^Oij + A 4 (0i^2)3 1 ivtm)^ + A 5(<Ai^2)3 2 (vtv2)3 2 ( A - 10 ) 

+ A 6 ((/>i^2)3 3 (??3r/2)3 3 + A 7 (0i02)3 5 (r?i%)3 5 g + A 8 (0i</> 2 )3 5 (vlm)3 5 A 

+ ?i (fyfah^ViVih! + / 2( ( / , j^)3 2? 3^i tr?2 )3 2? 3 + tiifa&jh^vtm)^ 



k 



k. 



i(x t x)3 1 (^i??2)3 1 + kitfwxh^ViWrn)^ 



fi) Hi) 

are needed. Let us briefly outline how the various couplings act: The couplings k\ and l\ 

renormalize Mi, fcg 1 -' splits masses of charged and neutral components, Ai and A2 mix neutral 

scalar and pseudoscalar components of the various fields. Hence, it also splits the masses of 

(i) (i) 

scalar and pseudoscalar of the lightest mass eigenstate, k±, l\ , Z3 mix the components of 



the various fj and adds flavor breaking effects. Since (^x\ J = such couplings do give 

contributions to mass terms and are not shown here. A3, . . . , As break Z4 and therefore mix 
components of 773 with components of 771,2- 

B Group Theory 

In this section, we give a short review of the relevant group theory of Qs x A4. We give the 
presentation of the group and a possible set of generators for all irreducible representations 
of the group. We summarize the most important Clebsch-Gordan coefficients for the quartet 
4 j and triplets 3j. See [32] for a more detailed description of the group theory of Qs >i A4. 
All Clebsch-Gordan coefficients can be obtained with the help of the Mathematica package 
Discrete, which has been published as part of [32]. 

B.l Mini-Review 

The semidirect product Qs x A4 we are using is defined by the relations 

SXS~ 1 =X, SYS~ 1 = Y- 1 , TXT~ 1 = YX, TYT' 1 = X . (B.l) 
between the generators of A4 

(S, T\S 2 = T 3 = (ST) 3 = 1} , (B.2) 

and Qs 

(X, Y|X 4 = 1, X 2 = Y 2 , Y~ l XY = X- 1 ) . (B.3) 
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Il 12 13 3 X 3 2 3 3 3 4 3 5 4i 4 2 4 3 



5 111 5*3 T3S3T3 T3S3T3 I3 T3S3T3 54 5*4 5*4 

T 1 (J Cl> 2 T3 T3 T3 T3 T3 T4 CJ 2 T4 WT4 

x 111 i 3 s 3 ^SsTs r 3 s 3 r 3 2 T$S 3 T 3 X A X A X A 

Table 2: Representations of Q$ xi A4 in the chosen basis. The first 4 representations are the unfaithful 
A& = (S 1 , T) representations to which the leptons are assigned (with p(X) = 1). Note that the representations 4j 
are double valued, i.e. p(Z(G) — X 2 ) = — 1, whereas the other representations are sing le valued {p{X 2 ) = 1). 
1-2 3 and 42 3 are complex, the other representations are real. 



Note that it is sufficient to use e.g. the generators X, S, T as Y = T 1 XT. The defining 
representation matrices for the representations we are using are given in Tab. 2 with 






f 





1 






) - 












1 - 




\ 


1 










10 
10 



(B.4) 



and S4 = 03 ® 0"i and X 4 = — ia 2 <8> 03. 



B.2 Clebsch-Gordan CoefRcients: Quartets 

The most important Clebsch-Gordan coefficients for the quartets a, b ~ 4 ^ are given by: 



(a f 6)] 



1 



a[bi + a l 2 b 2 + a 3 6 3 + a 4 6 4 



,tj 



(B.5) 



and the triplets: 



(a f 6)i 



(ot6)i 



(at6)i 



/ -0461 + a^ 2 - 0363 + a\b 4 



—0361 
a^i - 



- a\b2 + a} 63 + a\bi 
a\b 2 - a\b s + a\b 4 

( a\bi — a\b 2 — 0363 + 0464 
-a\bi + a\b 2 — 0363 + 0464 
-a\bi - a\b 2 + a\b 3 + a\b± 



( -alb 
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a^6 2 + <4&3 + <4&4 
a^i + a{&2 — O4&3 + a^4 



^ a 4 &i + a|&2 + "4^3 + a^4 
0361 + a\b 2 + a\b3 + a\b± 
y <4&i + a{f> 2 + 0463 + a\b 4 

f a\b\ — a\b 2 — a\b^ + a\b 4 
\ 



a\b\ + a\b 2 - a\bs + a} 2 b 4 



-a\b\ 



a\b 2 + a\bs + a\b 4 
(B.6) 



Note that 



(b a ) 3. is real for i = 2,3,4 and (0^6)3. 



-(6t ( 



for i = 1,5. 
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B.3 Clebsch-Gordan Coefficients: Triplets 



Furthermore, the most important Clebsch-Gordan coefficients for the 3-dimensional represen- 
tations 3j are described by 



(^)l. 



-j=. (a\b\ + a\b 2 + a\b 3 
f 



Jb) 



A,3 



a\b\ + u^a^bi + ^4^3 J ( a ^)l 

( 4&3 - a\b 2 

a\b x - a\b 3 | (a t& )s,3 
y a[b 2 - a\b\ 



a\bi + ua\b 2 + u?a\b^j (B.7) 

( a\b 3 + a\b 2 
a\bi + a\b 3 
\ a\b 2 + a\bi 



where [a\, a 2 , 03), (b±, b 2 , 63) ~ 3. Note that (a^a) A 3 is imaginary and (a^a) s ^ is real. Other 
important products are the product of a ~ 35 and b ~ 3^: 



(0^)3! 



/ 4^3 

4&i 



(a f 6)i 



of a ~ 3g and b ~ 32 

/ 4&a 



4&3 
y 4&1 



/ ^362 

4 63 
V 4 & i 



/ 4^2 
4^3 

\ a\bi 



(^6)33 



/ 4 fe 3 
4 6l 

V 4 ft 2 



/ 4 b 3 
a\b\ 
\ a\b 2 



and of a ~ 34 and b ~ 32 



/ 46; 



(at&). 



2 W 3 



3 01 



V 4&2 




*t6). 




(B.8) 



(B.9) 



(B.IO) 
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